Thermodynamics of ferromagnetic superconductors with spin-triplet electron pairing 
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We present a general thermodynamic theory that describes phases and phase transitions of fer- 
romagnetic superconductors with spin-triplet electron Cooper pairing. The theory is based on 
extended Ginzburg-Landau expansion in powers of superconducting and ferromagnetic order pa- 
rameters. We propose a simple form for the dependence of theory parameters on the pressure that 
allows correct theoretical outline of the temperature-pressure phase diagram for which at low tem- 
peratures a stable phase of coexistence of p-wave superconductivity and itinerant ferromagnetism 
appears. We demonstrate that the theory is in an agreement with the experimental data for some 
intermetallic compounds that are experimentally proven to be itinerant ferromagnetic exhibiting 
spin-triplet superconductivity. Some basic features of quantum phase transitions in such systems 
are explained and clarified. We propose to group the spin-triplet ferromagnetic superconductors in 
two different types of thermodynamic behavior, on the basis of quantitative criterion deduced from 
the present theory and the analysis of experimental data. 

PACS numbers: 74.20.De, 74.25. Dw, 64.70.Tg 



I. INTRODUCTION 

The spin-triplet or p-wave pairing allows parallel spin 
orientation of the fermion Cooper pairs in superfluid 3 He 
and unconventional superconductors For this rea- 
son the resulting unconventional superconductivity is ro- 
bust with respect to effects of external magnetic field 
and spontaneous ferromagnetic order, so it may coex- 
ist with the latter. This general argument implies that 
there could be metallic compounds and alloys, for which 
the coexistence of spin-triplet superconductivity and fer- 
romagnetism may be observed. 

Particularly, both superconductivity and itinerant fer- 
romagnetic orders can be created by the same band elec- 
trons in the metal, which means that spin-1 electron 
Cooper pairs participate in the formation of the itinerant 
ferromagnetic order. Moreover, under certain conditions 
the superconductivity is enhanced rather than depressed 
by the uniform ferromagnetic order that can generate it, 
even in cases when the superconductivity does not appear 
in a pure form as a net result of indirect electron-electron 
coupling. 

The coexistence of superconductivity and ferromag- 
netism as a result of collective behavior of /-band elec- 
trons has been found experimentally for some Uranium- 
based intermetallic compounds as, UGe2 
URhGe [1 0, 1, UCoGe @, 03, and Ulr 
low temperature (T ~ 1 K) all these compounds exhibit 
thermodynamically stable phase of coexistence of spin- 
triplet superconductivity and itinerant (/-band) electron 
ferromagnetism (in short, FS phase). In UGe2 and Ulr 
the FS phase appears at high pressure (P ~ 1 GPa) 
whereas in URhGe and UCoGe, the coexistence phase 
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persists up to ambient pressure (10 5 Pa= lbar). 

Experiments, carried out in ZrZn2 [131 ] , also indicated 
the appearance of FS phase at T < 1 K in a wide range of 
pressures (0 < P ~ 21 kbar). In Zr-based compounds the 
ferromagnetism and the p-wave superconductivity occur 
as a result of the collective behavior of the d-band elec- 
trons. Later experimental results [li], 03 had imposed 
the conclusion that bulk superconductivity is lacking in 
ZrZn 2 , but the occurrence of a surface FS phase at sur- 
faces with higher Zr content than that in ZrZn 2 has been 
reliably demonstrated. Thus the problem for the coex- 
istence of bulk superconductivity with ferromagnetism 
in ZrZn 2 is still unresolved. This raises the question 
whether the FS phase in ZrZn 2 should be studied by 
surface thermodynamics methods, or it can be investi- 
gated by considering that bulk and surface thermody- 
namic phenomena can be treated on the same footing. 
Taking in account the mentioned experimental results 
for ZrZnoand their interpretation by the experimental- 
ists m ei m we assume that the unified thermody- 
namic approach can be applied, although some specific 
properties of the FS phase in ZrZn 2 , if its surface nature 
is confirmed, will certainly need special study by surface 
thermodynamics. 

Here we will investigate the itinerant ferromagnetism 
and superconductivity of U- and Zr-based intermetal- 
lic compounds within the same general thermodynamic 
approach. Arguments supporting our point of view are 
given in several preceding investigations. We should men- 
tion that the spin-triplet superconductivity occurs not 
only in bulk materials, but also in quasi-two-dimensional 
(2D) systems - thin films and surfaces and quasi-lD wires 
(see, e.g., Refs. [16]). In ZrZn 2 both ferromagnetic and 
superconducting orders vanish at the same critical pres- 
sure P c , a fact implying that the respective order param- 
eter fields strongly depend on each other and should be 
studied on the same thermodynamic basis [T^. 
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The general thermodynamic treatment does not neces- 
sarily specify the system spatial dimensionality D: D = 3 
describes the bulk properties, and D = 2 - very thin films 
and mono- atomic layers. Within Landau theory of phase 
transitions (see, e.g., Ref. |l8|), the system dimension- 
ality can be distinguished by the values of the Landau 
parameters. Here we specify the values of these parame- 
ters from the experimental data for spin-triplet ferromag- 
netic superconductors. When the Landau parameters are 
obtained from microscopic theories, their values depend 
on the dimension D, only if the respective theory takes 
into account relevant fluctuation modes of order parame- 
ter fields, including long-scale fluctuation modes. We are 
not aware of well developed theories of this type which 
may figure out the complex behavior of the mentioned 
systems. Even in simple theories of band electron mag- 
netism, the Landau parameters are very complex func- 
tions of the density of states at the Fermi level and related 
microscopic parameters. Such complexity does not allow 
direct comparison between the results from microscopic 
theory and the experimental data. To make a progress in 
this situation we assume that the material parameters of 
our theory are loosely defined and may have values, cor- 
responding to various approximate microscopic theories, 
as mean-field approximation, spin-fluctuation theory, etc. 

For all compounds, cited above, the FS phase occurs 
only in the ferromagnetic phase domain of the T — P di- 
agram. Particularly at equilibrium, and for given P, the 
temperature Tp{P) of the normal-to- ferromagnetic phase 
(or N-FM) transition is never lower than the tempera- 
ture Tps(P) of the ferromagnetic-to-FS phase transition 
(FM-FS transition). This confirms the point of view that 
the superconductivity in these compounds is triggered by 
the spontaneous magnetization M , in analogy with the 
well-known triggering of the superfluid phase Ai in 3 He 
at mK temperatures by the external magnetic field H. 
Such "helium analogy" has been used in some theoretical 
studies (see, e.g., Ref. 0, H(|), where Ginzburg-Landau 
(GL) free energy terms, describing the FS phase were 
derived by symmetry group arguments. The non- unitary 
state, with a non-zero value of the Cooper pair magnetic 
moment, known from the theory of unconventional su- 
perconductors and superfluidity in 3 He [1], has been sug- 
gested firstly in Ref. , and later confirmed in other 
studies 0, [2(| ; recently, the same topic was comprehen- 
sively discussed in (2lj . 

For the spin-triplet ferromagnetic superconductors the 
trigger mechanism was recently examined in detail [22], 
[23( | . The system main properties are specified by the 
term in the GL expansion of form M\xp\ 2 , which repre- 
sents the interaction of M = {Mj;j — 1,2,3} with the 
complex superconducting vector field ip = {tpj}- Particu- 
larly, this term is responsible for the appearance of ip ^ 
for certain T and P values. A similar trigger mechanism 
is familiar in the context of improper ferroelectrics [2~jj . 

A crucial feature of these systems is the nonzero mag- 
netic moment of the spin-triplet Cooper pairs. As men- 
tioned above, the microscopic theory of magnetism and 



superconductivity in non-Fermi liquids of strongly in- 
teracting heavy electrons (/ and d band electrons) is 
either too complex or insufficiently developed to de- 
scribe the complicated behavior in itinerant ferromag- 
netic compounds. Several authors (see pjj H(| HH, l22l j) 
have explored the phenomenological description by a self- 
consistent mean field theory, and we will essentially use 
the thermodynamic results, in particular, the analysis 
in Refs. [22|, Mean-field microscopic theory of spin- 
mediated pairing leading to the mentioned non-unitary 
superconductivity state has been developed in Ref. [l7| 
that is in conformity with the phenomenological descrip- 
tion that we have done. 

In this paper, we present general thermodynamic treat- 
ment of systems with itinerant ferromagnetic order and 
superconductivity due to spin-triplet Cooper pairing of 
the same band electrons, which are responsible for the 
spontaneous magnetic moment. We outline their T — P 
phase diagrams and demonstrate two contrasting types of 
thermodynamic behavior. The present phenomenological 
approach includes both mean-field and spin-fluctuation 
theory (SFT), as the arguments in Ref. [25|. We propose 
a simple, yet comprehensive, modelling of P dependence 
of the free energy parameters, resulting in a very good 
compliance of our theoretical predictions for the shape 
the T — P phase diagrams with the experimental data 
(for some preliminary results, see Ref. [26j |. 

The theoretical analysis is done by the standard meth- 
ods of phase transition theory [l^]. Treatment of fluc- 
tuation effects and quantum correlations [H, [27| is not 
included in this study. But the parameters of the gener- 
alized GL free energy may be considered either in mean- 
field approximation as here, or as phcnomcnologically re- 
normalized parameters which are affected by additional 
physical phenomena, as for example, spin fluctuations. 

We demonstrate with the help of present theory that 
we can outline different possible topologies for the T — P 
phase diagram, depending on the values of Landau pa- 
rameters, derived from the existing experimental data. 
We show that for unconventional (spin-triplet) ferromag- 
netic superconductors (FSs) there exist two distinct types 
of behavior, which we denote as Zr-type (or, alternatively, 
type I) and U-type (or, type II). This classification of 
the FS, first mentioned in Ref. [26|], is based on the re- 
liable interrelationship between a quantitative criterion 
derived by us and the thermodynamic properties of the 
spin-triplet FSs (see Sec. III.C and Sec. IV. D). Our ap- 
proach can be also applied to URhGe, UCoGe, and Ulr. 
Our results shed light on the problems connected with 
the order of the quantum phase transitions at ultra-low 
and zero temperatures. They also raise the question for 
further experimental investigations of the detailed struc- 
ture of the phase diagrams in the high-P/low-T region. 

In Sec. II we present the GL free energy of ferromag- 
netic superconductors with spin-triplet Cooper pairing. 
Here the arguments for the proposed simple pressure de- 
pendence of the Landau parameters are explained and, 
as the comparison with experimental data demonstrates, 
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the chosen P-dependence gives quite satisfactory quanti- 
tative results. In more detailed studies one may use more 
complex pressure dependence of material parameters. In 
Sec. Ill we present our theoretical results for the vari- 
ous possible shapes of the phase diagram, multi-critical 
points, quantum phase transitions. The theory predicts 
all types of phase transition points and lines which are 
needed in the description of real itinerant ferromagnets 
with spin-triplet Cooper pairing. In Sec. IV we com- 
pare our theoretical results with the experimental data 
for some intermetallic compounds. The comparison gives 
very good quantitative agreement between experiment 
and theory for all known experimental examples, except 
URhGe, where a particular form of pressure dependence 
of the material parameters should be suggested. On the 



basis of this analysis we propose a possible classification 
of unconventional ferromagnetic superconductors, show- 
ing on the basis of quantitative criterion that they are 
of two types (Sec. IV. E). In Sec. V we summarize our 
findings and some unresolved problems which are beyond 
the scope of the present study. 



II. GENERALIZED GINZBURG-LANDAU FREE 
ENERGY 



Following Ref. |22|], the free energy per unit volume, 
F/V = f(ip, M), can be written in the form 



J 
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- a f M 2 + yM 4 + i 7o M • O x */>*) + SM 2 \ip\ 2 . (1) 
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The material parameters satisfy b s > 0, 6/ > 0; a s — 
a s (T-T s ), and a f = a f [T n - TJ(P)], where n = 1 gives 
the standard form of a/ , and n — 2 applies for SFT [25j 
and the Stoner-Wohlfarth model [28[ . The terms propor- 
tional to u s and v s describe, respectively, the anisotropy 
of the spin-triplet electron Cooper pairs and the crystal 
anisotropy. Next, 70 ~ J (with J > 0, the ferromag- 
netic exchange constant) and 6 > are parameters of 
the ip-M interaction terms. Previous studies [13] have 
shown that the anisotropy represented by the u s and v s 
terms in Eq. ([1]), slightly perturbs the size and shape of 
the stability domains of the phases, while similar effects 
can be achieved by varying the b s factor in the 6 s |i/'| 4 
term. For such reasons, in the present analysis we ignore 
the anisotropy terms, setting u s = v s — 0, and consider 
b s = b > as an effective parameter. Then, without loss 
of generality, we may choose the magnetization vector to 
have the form M = (0, 0, M). 

According to the microscopic theory of band mag- 
netism and superconductivity, the macroscopic material 
parameters in Eq. ([1} depend in a very complex way on 
the density of states at the Fermi level and related mi- 
croscopic quantities [2^ | . That is why we can hardly use 
the microscopic characteristics of these complex metallic 
compounds in order to elucidate their thermodynamic 
properties, in particular, in outlining their phase dia- 
grams in some details. However, some simple microscopic 
models give useful results, for example, the zero temper- 
ature Stoner-type model employed in Ref. [30j . 

We redefine for convenience the free energy, Eq. (fT]), 
in a dimensionless form by / = //(6/Mq), where Mq = 
[afTJ^/bf] 1 / 2 > is the value of the magnetization M 
corresponding to the pure magnetic subsystem (i/> = 0) 



at T = P = and T/ = T/(Q). Order parameters 
assume the scaling m = M/Mq and ip = ■0/[(6y/6) 1 / 4 Mo] 
and in result, the free energy becomes 



/ = rcj) 2 + — +tm 2 -\ h27TO</> 1 2 sin0 + 7i?7r</r, (2) 



where 



1) 

1O1 



h = \<Pjl (j = ijM), 4> = \<f\, and 9 = yu 2 
is the phase angle between the complex ipi — <p\e 
and ip 2 — 02e 02 - Note, that the third phase angle, 
6*3, corresponding to the third complex field component 



V3 



does not enter explicitly in the free energy 
density /, given by Eq. @, which is a natural result of 
the continuous space degeneration. The dimensionless 
parameters t, r, 7 and 71 in Eq. ^ are given by 



t = T n - Tf(P) 
where k = a s by 2 /ctf^^T 

71 = s/{bb s y/ 2 



= k(T-T s ), 



(3) 



/V 1 , 7 = l /[a f T?S y \ _a*d 
The reduced temperatures are: T = 
T/T /0 , T f (P) - T f (P)/T f0 , and f s (P) =T S (P)/T f0 . 

Here we will outline the simple assumptions for the P- 
dependence of the i, r, 7, and 71 parameters in Eq. 
Particularly, we assume that only Tf has a significant 
P-dependence, described by 



T f (P) = (l-P) 



1/,. 



(4) 



where P = P/Po and Po is a characteristic pressure de- 
duced later. In ZrZn2 and UGe2 the Po values are very 
close to the critical pressure P c , at which both the fer- 
romagnetic and superconducting orders vanish, but in 
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other systems this is not necessarily the case. As we 
will discuss, the nonlinearity (n — 2) of P/(P) in ZrZn 2 
and UGe2 is relevant at relatively high P, at which the 
N-FM transition temperature Tp(P) may not coincide 
with T f (P). The form J3]) of the model function f f (P) 
is consistent with preceding experimental and theoreti- 
cal investigations of the N-FM phase transition in ZrZn 2 
and UGe2 (see, e.g., Refs. 0, Hfl Here we consider 

only non-negative values of the pressure P; for effects at 
P < 0, see, e.g., Ref. [H|. 

The model function (UJ) is defined for P < P , in par- 
ticular, for the case of n > 1, but we should have in mind 
that, in fact, the thermodynamic analysis of Eq. (J2]) in- 
cludes the parameter t rather than T/(P). This param- 
eter is given by 



t(T, P) = T n - 1 + P, 



(5) 



and is well defined for any P. This allows to consider 
pressures P > P within the free energy ©. 

The model function P/(P) can be naturally general- 
ized to f/(P) = (1 - P 3 ) 1 /" but the present needs of 
interpretation of experimental data do not require this; 
hereafter we use Eq. © which corresponds to j3 = 1 and 
a — n). Besides, other analytical forms of T/(P) can also 
be tested in the free energy ([2|), in particular, expansion 
in powers of P, or, alternatively, in (1 — P) which sat- 
isfy the conditions f f (0) = 1 and f/(l) = 0. Note, that 
in URhGe the slope of T F (P) ~ T/(P) is positive from 
P = up to high pressures [f| and for this compound 
the form © of P/(P) is inconvenient. Here we apply the 
simplest form for the P-dependence, namely, Eqs. ^ - 
©■ 



III. PHASE TRANSITIONS AND 
MULTICRITICAL POINTS 

A. Phase of coexistence of ferromagnetism and 
superconductivity 

The simplified model © is capable of describing the 
main thermodynamic properties of spin-triplet ferromag- 
netic superconductors. For r > 0, i.e., T > T s , there 
are three stable phases: (i) the normal (N-) phase, given 
by <j) = m — (stability conditions: t > 0, r > 0); 
(ii) the pure ferromagnetic phase (FM phase), given by 
to = (— i) 1 / 2 > 0, 4> — 0, which exists for t < and is 
stable, provided r > and r > (71 1 + 7|£| 1 / 2 ), and (iii) 
the phase of coexistence of ferromagnetic order and su- 
perconductivity (FS phase), given by sin# = =pl, fa = 0, 
fa = fa = 4>/ \/2, where 



<j) 2 = k(T s - f) ± 7 to - 71 to 2 > 0. (6) 
The magnetization to satisfies the equation 



C3TO 3 ± C2TO 2 + ClTO ± Co = 

with coefficients cq = 7«(P — T s ), 



ci 



f n + K-/ 1 (f s -f)+P-l 



(-2 



3771, c 3 = 2(l-7 2 ; 



(7) 



(8) 



(9) 



The FS phase has two thermodynamically equivalent 
phase domains that can be distinguished by the upper 
and lower signs (±) of some terms in Eqs. (|6|) and jTJ. 
The upper sign describes the domain (labelled bellow 
again by FS), with to > 0, sm9 = —1, whereas the lower 
sign describes the conjunct domain FS*, where to < 
and sin^ = f ; for details, see, Ref. J2^|. Here we consider 
one of the two thermodynamically equivalent phase do- 
mains, namely, the domain FS, which is stable for m > 
(FS* is stable for m < 0). In the absence of external 
symmetry breaking field, such "one-domain approxima- 
tion" correctly presents the main thermodynamic prop- 
erties described by the model {l}. The FS phase domain, 
given by upper sign Eqs.© and (J7J), satisfies the stability 
conditions: jM > and 



K(f s - f ) ± 7TO - 2 7 l 771 2 > 0, 



(10) 



3(l-7 2 )TO 2 + 377iTO + r n -l + P+ K7l (T s -T)-^- > 0. 

(11) 

These results are valid for TAP) > T S (P), which ex- 
cludes any pure superconducting phase (if) ^ 0, to = 0) 
in accordance with the available experimental data. 

For r < 0, and t > the model ©-© exhibits a 
stable pure superconducting phase (fa = fa = m = 0, 
2 = — r) 22]. This phase may occur in the temper- 
ature domain Tf(P) < T < T s . For systems, where 
T/(0) 3> T s , this is a domain of pressure in a very close 
vicinity of Pq ~ P c , where Tp(P) ~ Tf(P) decreases to 
values lower than T s . Such a situation is described by the 
model ([2]) only if T s > 0. This case is interesting from 
the experimental point of view only when T s > is big 
enough to enter in the scope of experimentally measur- 
able temperatures. Up to date a pure superconducting 
phase has not been observed within the accuracy of ex- 
periments on the studied metallic compounds. For this 
reason, in the reminder of this paper we will assume that 
the critical temperature T s of the generic superconduct- 
ing phase transition is either non-positive (T s < 0), or 
has a small positive value which can be neglected in the 
analysis of the available experimental data. 

The negative values of the critical temperature T s of 
the generic superconducting phase transition are gen- 
erally possible and produce a variety of phase diagram 
topologies (Sec. III.B). The value of T s depends on the 
strength of the interaction, responsible for the formation 
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TABLE I: Theoretical results for the location [(T , P)- reduced coordinates] of the tricritical points A = (Ta,Pa) and B = 
(Tb, Pb), the critical-end point C = (Tq, Pa), an d the point of temperature maximum, max = {T m , P m ) on the curve Tps(P) of 
the FM-FS phase transitions of first and second order (for details, see Sec. III.B). The first column shows Tn = T(A,B,c,m)' The 
second column stands for ijv = i(A,s,c,m) • The reduced pressure values P(A,B,c,m) of points A, B, C, and max are denoted by 
Pjv(n): n — 1 stands for the linear dependence Tf(P), n = 2 - for the nonlinear dependence of Tf(P) and t(T), corresponding 
to SFT. 



N 


f N 




P N (n) 


A 
B 
C 

max 


T a 

f s + 7 2 (2 + 7 i)/4 K (l+ 7l ) 2 

T s +7 2 /4k(1+7i) 

T s +7 2 /4k7i 


7 2 /2 

-7 2 /4(l + 7i) 2 


- 7 2 /4 7l 2 


1 - T s " +7 2 /2 

i-rS-7 2 /4(i + 7i) 2 
1-fS 

1 - T" - 7 2 /47? 



of the spin-triplet Cooper pairs, which means that the 
P-dependence of T s specifies the sensitivity of electron 
couplings to the crystal lattice properties. This is an ef- 
fect which might be included in our theoretical scheme by 
introducing some convenient pressure dependence of T s . 
To do so we need information either from experimental 
data or from a comprehensive microscopic theory. 

Usually, T s < is interpreted as a lack of any super- 
conductivity, but here the same non-positive values of T s 
are effectively enhanced to positive values by the inter- 
action parameter 7 which triggers the superconductiv- 
ity up to superconducting phase transition temperatures 
Tfs(P) > 0. This is readily seen from Table 1, where 
we present the reduced characteristic temperatures on 
the FM-FS phase transition line T F s(P), calculated from 
the present theory, namely, T m - the maximum of the 
curve Tps(P) (if available; see Sec. III.B), the tempera- 
tures Ta and Tb, corresponding to the tricritical points 
A = (f A ,P A ), B = (f B ,P B ), and the temperature f c 
corresponding to the critical-end point C = (Tc,Pc)- 
The theoretical derivation of the dependence of the multi- 
critical temperatures Ta, Tb, and Tq on 7, 71, k, and T s , 
as well as the dependence of T m on the same model pa- 
rameters is outlined in Sec. III.B. All these temperatures 
as well as the whole phase transition line Tps(T) are 
considerably boosted above T s owing to positive terms 
of order 7 2 . If T s < 0, the superconductivity appears, 
provided T m > 0, i.e., when 7 2 /4^7i > \T S \ (see Table 
1)- 



B. Temperature-pressure phase diagram 

Although the structure of the FS phase is complicated, 
some of the results can be obtained in analytical form. 
A more detailed outline of the phase domains, for exam- 
ple, in T — P phase diagram, can be done by using suit- 
able values of material parameters in the free energy 
Pq, Tfo, T s , k, 7, and 71. Here we present some of the 
analytical results for the phase transition lines and the 
multi-critical points. Typical shapes of phase diagrams 
derived directly from Eq. ^ are given in Figs. 1-6. Fig. 1 
shows the phase diagram calculated from Eq. @ for pa- 



rameters, corresponding to the experimental data [13l ] for 
ZrZn2. Figs. 2-3 show the low-temperature and the high- 
pressure parts of the same phase diagram (see Sec. IV 
for details). Figs. 4-6 show the phase diagram calculated 
for the experimental data 0, Q of UGe2 (see Sec. IV). 
In ZrZn2, UGe2, as well as in UCoGe and Ulr, critical 
pressure P c exists, where both superconductivity and fer- 
romagnetic orders vanish. 

As in experiments, we find out from our calculation 
that in the vicinity of To ~ T c , the FM-FS phase tran- 
sition is of first order, denoted by the solid line BC in 
Figs. 2,3,5,6. At lower pressure the same phase tran- 
sition is of second order, shown by the dotted lines in 
the same figures. The second order phase transition line 
Tfs(P) separating the FM and FS phases is given by the 
solution of the equation 

Tfs(P) =f s + fi*Fs(P) + i[-t FS (P)] 1/2 , (12) 

where t FS (P) =t(T FS ,P) < 0, 7 = j/k, 71 = ji/n, 
and < T < Pb; Pb is the pressure corresponding 
to the multi-critical point B, where the first-order line 
Tps(P) terminates, as clearly seen in Figs. 3 and 6. 
Eq. p2[) strictlv coincides with the stability condition for 
the FM-phase with respect to the appearance of FS- 
phase. 

Additional information for the shape of the FM-FS 
phase transition line can be obtained by the derivative 
p = df FS (P)/dP, namely, 

Ps+ii-m-tps) 1 ' 2 

where p s = df s (P)/dP. Note, that Eq. is obtained 
from Eqs. © and ^HJ). 

The shape of the line T F $(P) can vary depending on 
the theory parameters (see, e.g., Figs. 2, 5). For certain 
ratios of 7, 71, and values of p s , the curve T F s(P) ex- 
hibits a maximum T m = T F s(P m ), given by the solution 
of p(p s ,T m , P m ) = 0. This maximum is clearly seen in 
Figs. 5, 6. To locate the maximum we need to know 
p s . We have already assumed that T s does not depend 
on T, as explained above, which from the physical point 
of view means that the function T S (P) is flat enough to 
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10 15 

(kbar) 



FIG. 1: T - P diagram of ZrZn 2 calculated for T a = 0, T f0 = 
28.5 K, P = 21 kbar, k = 10, 7 = 2fi « 0.2, and n = 1. The 
dotted line represents the FM-FS transition and the dashed 
line stands for the second order N-FM transition. The dotted 
line has a zero slope at P = 0. The low-temperature and 
high-pressure domains of the FS phase are seen more clearly 
in the following Figs. 2, 3. 



allow the approximation T s rs without a substantial 
error in the results. From our choice of P-dependence of 
the free energy ©parameters, it follows that p s = 0. 
Setting p s = p = in Eq. (fTS"]) , we obtain 
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t(T m , Pm) 4^y 2 ' 



(14) 



namely, the value t m (T, P) = t(T m , P m ) at the maximum 
Tm{Pm) of the curve Tps(P). Substituting t m back in 
Eq. (|12p we obtain T m , and with its help we also obtain 
the pressure P m , both given in Table 1, respectively. 

We want to draw the attention to a particular feature 
of the present theory that the coordinates T m and P m of 
the maximum (point "max") on the curve Tps(P), and 
also the results of various calculations with the help of 
Eqs. (JT3J) and (TP3")) are expressed in terms of the reduced 
interaction parameters 7 and 71. Thus, using certain 
experimental data for T m , P m , as well as Eqs. (|12p and 
(|13p for Tps, T s and the derivative p at particular values 
of the pressure P, 7 and 71 can be calculated without any 
additional information, for example, for the parameter 
k. This property of the model © is very useful in the 
practical work with the experimental data. 

The conditions for existence of a maximum on the 
curve Tps{P) can be determined by requiring P m > 
and T m > and using the respective formulae for these 
quantities, shown in Table 1. This max always occurs 
in systems where Tp$(0) < and the low-pressure part 
of the curve Tps(P) terminates at T = for some non- 
negative critical pressure Pq c (see Sec. III.C). But the 
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FIG. 2: Details of Fig. 1 with expanded temperature scale. 
The points A, B, C are located in the high-pressure part 
(P ~ P c ~ 21 kbar). The max point is at P ~ kbar. 
The FS phase domain is shaded. The dotted line shows the 
second order FM-FS phase transition with P m m 0. The solid 
straight line BC shows the hst-order FM-FS transition for 
P > Pb- The quite flat solid line AC shows the first order 
N-FS transition (the lines BC and AC are more clearly seen 
in Fig. 3. The dashed line stands for the second order N-FM 
transition. 



maximum may occur also for some sets of material pa- 
rameters, when Tpg(Q) > (see Fig. 2, where P m = 0). 
All these shapes of the line Tps(P) are described by the 
model ((2|). Irrespective of the particular shape, the curve 
Tps(P) given by Eq. ([12")) always terminates at the tricrit- 
ical point (labelled B), with coordinates (Pb,Tb) (see, 
e.g., Figs. 3, 6). 

At pressure P > Pb the FM-FS phase transition is of 
first order up to the critical-end point C. For Pb < P < 
Pc the FM-FS phase transition is given by the straight 
line BC; see, Figs. 3, 6. All three phase transition lines 
- N-FM, N-FS, and FM-FS - terminate at point C. For 
P > Pc the FM-FS phase transition occurs on a rather 
flat, smooth line of equilibrium transition of first order up 
to a second tricritical point A with Pa ~ Po and Ta ~ 0. 
Finally, the third transition line terminating at the point 
C describes the second order phase transition N-FM. The 
reduced temperatures Tjq and pressures Pjv, N = (A, B, 
C, max) at the three multi-critical points (A, B, C) and 
the maximum X' m (P m ) are given in Table 1. Note that, 
for any set of material parameters, Ta < Tc < Tb < T m 
and P m < P B < P c < Pa- 

There are other types of phase diagrams, resulting from 
the model For negative values of the generic super- 
conducting temperature T s several other topologies of the 
T — P diagram can be outlined. The results for the multi- 
critical points, presented in Table 1 show that when T s 
lowers below T = 0, also Tc decrease, firstly to zero, 
and then to negative values. When Tc = the direct 
N-FS phase transition of first order disappears and the 
point C becomes a very special zero-temperature multi- 
critical point. As seen from Table 1, this happens for 
T s = -7 2 T/(0)/4k(1 + 71). The further decrease of T s 
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FIG. 3: High-pressure part of the phase diagram of ZrZri2, 
shown in Fig. 1. The thick solid lines AC and BC show the 
first-order transitions N-FS, and FM-FS, respectively. Other 
notations are explained in Figs. 1, 2. 



causes the point C to fall below the zero temperature and 
then the zero-temperature phase transition of first order 
near P c splits into two zero-temperature phase transi- 
tions: a second order N-FM transition and a first order 
FM-FS transition, provided Tg still remains positive. 

At lower T s also the point B falls below T = and 
the FM-FS phase transition becomes entirely of second 
order. For very extreme negative values of T s , a very 
large pressure interval below P c may occur where the FM 
phase is stable up to T = 0. Then the line Tps(P) will 
exist only for relatively small pressure values (P <C P c ). 
This shape of the stability domain of the FS phase is also 
possible in real systems. 



C. Quantum phase transitions 

The effects of quantum correlations on the phase 
transition properties at ultra-low and zero temperature, 
known as quantum phase transitions [27j , can also be de- 
scribed by the free energy (Q]), as mentioned above. The 
time-dependent quantum fluctuations, which describe 
the intrinsic quantum dynamics of spin-triplet ferromag- 
netic superconductors, are not included in the model |T]), 
but some basic properties of the quantum phase transi- 
tions can be outlined within the classical limit, namely, 
by the free energy model (HJ) — - 

Generally, both classical (thermal) and quantum fluc- 
tuations are investigated by the method of the renormal- 
ization group (RG) [l8|, which is specially intended to 
treat the generalized action of system, where the order 
parameter fields fluctuate in time t and space x [H, H3 • 
These effects, which are beyond the scope of the paper, 
lead either to precise treatment of the narrow critical 
region in a very close vicinity of second order phase tran- 
sition lines, or to a fluctuation-driven change of the phase 



transition order. But, the thermal fluctuation and quan- 
tum correlation effects on the thermodynamics of a given 
system can be unambiguously estimated only after the 
results from the counterpart simpler theory, where these 
effects are not present, are known and, hence, the dis- 
tinction in the thermodynamic properties predicted by 
the respective variants of the theory can be established. 
Here we show that the basic low and ultra-low temper- 
ature properties of the spin-triplet FSs, as given by the 
preceding experiments, are derived from the model (1) 
without any account of fluctuation phenomena and quan- 
tum correlations. The latter might be of use in a more 
detailed consideration of the close vicinity of quantum 
critical points in the phase diagrams of spin-triplet FSs. 
Here we show that the theory predicts quantum critical 
phenomena only for quite particular physical conditions 
whereas the low- and zero-temperature phase transitions 
of the first order are favored by both symmetry argu- 
ments and detailed thermodynamic analysis. 

There is a number of experimental [t| [33| and theo- 
retical [13, [H, [H| studies on the problem of quantum 
phase transitions in unconventional ferromagnetic super- 
conductors, including the mentioned intermetallic com- 
pounds. Some of them are based on different theoretical 
schemes, and do not refer to the model (JTJ). Other, for ex- 
ample, those in Ref. UK reported results about the ther- 
mal and quantum fluctuations, described by the model 
© before the comprehensive knowledge of the results 
from the basic treatment, given by the present investi- 
gation. In such cases one could not be sure about the 
correct interpretation of the results from the RG and the 
possibilities for their application to particular zero tem- 
perature phase transitions. Here we present basic results 
for the zero-temperature phase transitions, described by 
the model ©. 

The RG investigation [3J| has demonstrated up to two 
loop order of the theory that the thermal fluctuations 
of the order parameter fields re-scale the model in 
way which corresponds to first order phase transitions 
in magnetically anisotropic systems. This result is im- 
portant for the metallic compounds, we consider here, 
because in all of them magnetic anisotropy is present. 
The uniaxial magnetic anisotropy in ZrZn2 is much more 
weak than in UGe2 but cannot be neglected when fluc- 
tuation effects are accounted for. Owing to the particu- 
lar symmetry of the model for the case of magnetic 
isotropy (Heisenberg symmetry), the RG study reveals an 
entirely new class of (classical) critical behavior. Besides, 
the different spatial dimensions of the superconducting 
and magnetic quantum fluctuations imply a lack of stable 
quantum critical behavior even when the system is com- 
pletely magnetically isotropic. The pointed arguments 
and preceding results lead to the reliable conclusion that 
the phase transitions, which have already been proven to 
be first order in the lowest order approximation, where 
thermal and quantum fluctuations are neglected, will not 
undergo a fluctuation-driven change of the phase transi- 
tion order from first to second. Such picture is described 
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FIG. 4: T - P diagram of UGe 2 calculated taking T 3 = 0, 
T f0 = 52 K, P = 1-6 GPa, k = 4, 7 = 0.0984, 71 = 0.1678, 
and n = 1. The dotted line represents the FM-FS transition 
and the dashed line stands for the N-FM transition. The low- 
temperature and high-pressure domains of the FS phase are 
seen more clearly in the following Figs. 5, 6. 



below, in Sec. IV, and it corresponds to the behavior of 
real compounds. 

Our results definitely show that the quantum phase 
transition near P c is of first order. This is valid for the 
whole N-FS phase transition below the critical-end point 
C, as well as the straight line BC. The simultaneous ef- 
fect of thermal and quantum fluctuations do not change 
the order of the N-FS transition, and it is unlikely to sup- 
pose that thermal fluctuations of the superconductivity 
field i/> can ensure a fluctuation-driven change of the or- 
der of the FM-FS transition along the line BC. Usually, 
the fluctuations of ip in low temperature superconduc- 
tors are small and slightly influence the phase transition 
in the very narrow critical region in the vicinity of phase 
transition point. This effect is very weak and can hardly 
be observed in any experiment on low-temperature su- 
perconductors. Besides, the fluctuations of the magnetic 
induction B always tend to a fluctuation-induced first or- 
der phase transition rather than to the opposite effect - 
the generation of magnetic fluctuations with infinite cor- 
relation length at the equilibrium phase transition point 
and, hence, a second order phase transition [lj|, l36i |. 
Thus we can reliably conclude that the first order phase 
transitions at low-temperatures, represented by the lines 
BC and AC in vicinity of P c do not change their order 
as a result of thermal and quantum fluctuations. 

Quantum critical behavior for continuous phase transi- 
tions in spin-triplet ferromagnetic superconductors with 
magnetic anisotropy can therefore be observed at other 
zero-temperature transitions, which may occur in these 
systems far from the critical pressure P c . This is possible 
when Tps(0) = and the Tps(P) curve terminates at 
T = at one or two quantum (zero-temperature) crit- 
ical points: Pq c < P m - "lower critical pressure" , and 



FIG. 5: Low-temperature part of the T — P phase diagram 
of UGe2, shown in Fig. 4. The points A, B, C are located in 
the high-pressure part (P ~ P c ~ 1.6 GPa). The FS phase 
domain is shaded. The thick solid lines AC and BC show the 
first-order transitions N-FS, and FM-FS, respectively. Other 
notations are explained in Figs. 1,2. 



Pq c > P m - "upper critical pressure." To obtain these 
critical pressures one should solve Eq. (TT2|) with respect 
to P, provided T FS {P) = 0, T m > and P m > 0, 
namely, when the continuous function Tps(P) exhibits 
a maximum. The critical pressure Pq c is bounded in the 
relatively narrow interval (P m ,Pe) and can appear for 
some special sets of material parameters (r, £, 7,71). In 
particular, as our calculations show, Pq c do not exist for 
T s > 0. 

The analytical calculation of the critical pressures Pq c 
and Pq c for the general case of T s ^ leads to very com- 
plex conditions for the appearance of second critical field 
Pq c . The correct treatment for T s ^ can be performed 
within the two-domain picture of the phase FS. The com- 
plete study of this case is beyond our aims, but here we 
will illustrate our arguments by investigation of the con- 
ditions, under which the critical pressure P oc occurs in 
systems with T s w 0. Moreover, we will present the gen- 
eral result for Pq c > and Pq c > in systems where 
T s ± 0. 

Setting T FS (P 0c ) = in Eq. ([H]) we obtain the follow- 
ing quadratic equation 



7im 0c - jm 0c - T s = 
for the reduced magnetization 

= H(o,p oc )] 1 / 2 = (i-p 0c ) 1/2 



(15) 



(16) 



rriQ C 



and, hence, for P 0c . For T s ^ the Eqs. ([15]) and HU) 
have two solutions with respect to Pq c . For some sets 
of material parameters these solutions satisfy the phys- 
ical requirements for Pq c and Pq c and can be identified 
with the critical pressures. The conditions for existence 
of Pbc and Pq c can be obtained either by analytical cal- 
culations, or by numerical analysis for particular values 
of the material parameters. 
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For T s = 0, the trivial solution Pq c = 1 corresponds to 
Poc = Po > Pb and, hence, does not satisfy the physical 
requirements. The second solution 



Poc = 1 - 



if 



is positive for 



2i>i 

7 



(17) 



(18) 



and, as shown below, it gives the location of the quantum 
critical point (T = 0,Pq c < P m )- At this quantum criti- 
cal point, the equilibrium magnetization moc is described 
by mo c = 7/71 and is twice bigger that the magnetization 
m m — 7/271 [HJ at the maximum of the curve Tps(P). 

By taking P m from Table 1, we can show that the 
solution (fT7| satisfies the condition P 0c < P m for n = 1, 
if 



71 < 3k 

and for n = 2 (SFT case), when 
7 < 2V3k. 



(19) 



(20) 



Finally, we determine the conditions under which the 
maximum T m of the curve Tpg(P) occurs at non- negative 
pressures. For n — 1, we obtain that P m > for n = 1, 
if 



7i 



> 

7 - 2 



whereas for n = 2, the condition is 



(21) 



li>i(l + Jl 
7 ~ 2 V 4k 2 



1/2 



(22) 



Obviously, the conditions (|18| - ([22[) are compatible 
with one another. The condition (|21[) is weaker than the 
condition (fl"8jl. provided the inequality (fl9|) is satisfied. 
The same is valid for the condition (|2"2")l . if the inequality 
([20"]) is valid. In Sec. IV we will show that these theoret- 
ical predictions are confirmed by the experimental data. 

Doing in the same way the analysis of Eq. (fT2"]l some 
results may easily be obtained for T s ^ 0. In this more 
general case the Eq. (fl2|) has two nontrivial solutions that 
yield two possible values of the critical pressure. 



(±) 



47? 



i± 1 



4T s K7i 
7 2 



1/2- 



(23) 



The relation -Fbe(-) > Pqc(+) is always true. Therefore, 
to have both P(w±) ^ 0: it is enough to require that 
Pqc(+) > 0. Having in mind that for the phase diagram 
shape, we study, T m > 0, and according to the result for 
T m in Table 1, this leads to inequality T s > — 7 2 /4«7i. 
So, we obtain that Pq c (+) > will exist, if 



£>1 
7 



kT s 

7 



(24) 



which generalizes the condition (fTS)) . 

Now we can identify the pressure Po c (+) with the lower 
critical pressure Pq c , and Poc(-) with the upper critical 
pressure Pq c . Therefore, for wide variations of the ma- 
terial parameters, the theory ([1]) describes a quantum 
critical point P oc that exists, provided the condition ([24]) 
is satisfied. The quantum critical point (T = 0, P 0c ) ex- 
ists in UGe2 and, perhaps, in other p-wave ferromagnetic 
superconductors, for example, in Ulr. 

Our results predict the appearance of second critical 
pressure - the upper critical pressure P' oc that exists un- 
der more restricted conditions and, hence, can be ob- 
served for those systems, where T s < 0. As mentioned 
in Sec. III.B, for very extreme negative values of T Sl 
when Tb < 0, the upper critical pressure Pq c > oc- 
curs, whereas the lower critical pressure Pq c > does 
not appear. But especially this situation should be in- 
vestigated in a different way, namely, one should keep 
Tfs{0) different from zero in Eq. (fl"2")) . and consider a 
form of the FS phase domain in which the curve Tps (P) 
terminates at T — for Pq c > 0, irrespective of whether 
the maximum T m exists or not. In such geometry of the 
FS phase domain, the maximum T^P m ) may exist only 
in quite unusual cases, if it exists at all. 

The quantum and thermal fluctuation phenomena in 
the vicinities of the two critical pressures Pq c and Pq c 
need a non-standard RG treatment because they are re- 
lated with the fluctuation behavior of the superconduct- 
ing field far below the ferromagnetic phase transitions, 
where the magnetization M does not undergo significant 
fluctuations and can be considered uniform. The pres- 
ence of uniform magnetization produces couplings of M 
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and tp which are not present in previous RG studies and 
need a special analysis. 



IV. APPLICATION TO METALLIC 
COMPOUNDS 

A. Theoretical outline of the phase diagram 

In order to apply the above displayed theoretic calcu- 
lations, following from the free energy @, for the outline 
of T — P diagram of any material, we need information 
about the values of Po, T/o, T s , k 7, and 71. The temper- 
ature T/o can be obtained directly from the experimen- 
tal phase diagrams. The pressure Po is either identical 
or very close to the critical pressure P c , for which the 
N-FM phase transition line terminates at T ~ 0. The 
temperature T s of the generic superconducting transition 
is not available from the experiments because, as men- 
tioned above pure superconducting phase not coexisting 
with ferromagnetism has been not observed. This can 
be considered as an indication that T s is small and does 
not produce a measurable effect. So the generic super- 
conducting temperature will be estimated on the basis of 
following general arguments. For T/(P) > T s we must 
have T S (P) = at P > P c , where X)(P) < and for 
< P < Po, T s < Tq- Therefore for materials where 
Tc is too small to be observed experimentally, T s can be 
ignored. 

As far as the shape of FM-FS transition line is well 
described by Eq. (p~2|) we will make use of additional 
data from available experimental phase diagram for fer- 
romagnetic superconductors. For example, in ZrZn 2 
these are the observed values of Tps(0) and the slope 
po = [dT FS (P)/dP} = (T f o/Po)po at P = 0. For UGe 2 , 
where a maximum is observed on the phase transition 
line, we can use the experimental values for T m , P m , and 
also for Po c . The interaction parameters 7 and 71 are 
derived using Eq. (fl~2|) and the expressions for T m , P m , 
and po, see Table 1. The parameter k is chosen by fitting 
the expression for the critical-end point Tc- 

B. Diagram of ZrZri2 

Experiments for ZrZn 2 [l3| give the following values: 
T/o = 28.5 K, T FS (0) = 0.29 K, P - P c = 21 kbar. The 
curve Tp(P) ~ Tf(P) is almost a straight line, which 
directly indicates that n = 1 is adequate in this case 
for the description of the P-dependence. The slope for 
Tfs(P) at P = is estimated from the condition that 
its magnitude should not exceed T/o/P c ~ 0.014 as we 
have assumed that the transition line is straight one, so 
in result we have —0.014 < p < 0. This ignores the pres- 
ence of a maximum. The available experimental data for 
ZrZn2, do not give clear indications whether a maximum 
at (T m ,P m ) exists. If there is a maximum at P = this 
will mean po — 0, and a maximum with T m ~ Tps(0) 



and P m <C Po will give < po < 0.005. We may choose 
that po = 0, then 7 « 0.02 and 71 » 0.01, similar values 
for the parameters are obtained for any |po| < 0.003. 

The multi-critical points A and C cannot be distin- 
guished experimentally as the experimental accuracy fl3T ] 
is less than ~ 25 mK in the high-pressure domain (P ~ 
20 — 21 kbar). Having this mind we can propose as ini- 
tial assumption that Tc ~ 10 mK, which corresponds to 
k ~ 10, see Table 1. We have calculated the T — P dia- 
gram using po = and po = 0.003, as well as the other 
parameters, discussed above. The differences, obtained 
in these two cases, are negligible. 

In Fig. 1 we show the phase diagram of ZrZn 2 cal- 
culated directly from the free energy ^ for n = 1, the 
above mentioned values of T s , Po, T/o, K, and values of 
7 « 0.2 and 71 « 0.1 which ensure p « 0. The ob- 
tained coordinates of characteristic points are the follow- 
ing: P A ~ P c = 21.42 kbar, T A = T F (P C ) = T FS (P C ) = 
K; P B = 20.79 kbar, T B = 0.0495 K, and P c = 20.98 
kbar, T c = 0.0259 K . 

In Fig. 2, the low-T region is seen in more detail with 
A, B, C points, respectively. The dotted Tps(P) curve, 
represents the phase transition line FM-FS of second or- 
der and has a maximum T m = 0.290 K at P = 0.18 kbar, 
which is slightly above T FS (0) = 0.285 K. The straight 
solid line BC in Fig. 2 shows the first order FM-FS phase 
transition which occurs for P B < P < Pc- The solid 
AC line shows the first order N-FS phase transition and 
the dashed line stands for the N-FM phase transition of 
second order. 

Although the expanded temperature scale in Fig. 2, 
the difference [T m — T F g(0)] = 5 mK can hardly be seen. 
The calculations, both analytical and numerical, lead to 
the conclusion that if we want to locate the point max 
exactly at P = we have to take values of 7 and 71 
with accuracy up to 10" 4 . This may be considered as 
an indication that the max for parameters correspond- 
ing to ZrZn 2 is very sensitive to very small changes of 
parameters 7 and 71 around the values 0.2 and 0.1, re- 
spectively. Our initial idea was to present a diagram 
with T m = T FS {0) = 0.29 K and po = 0, with the point 
of maximum being exactly at P = 0, but the final phase 
diagram slightly departs from this picture because of the 
mentioned sensitivity of the result on the values of the 
interaction parameters 7 and 71. 

We have also calculated the theoretical phase diagram 
of ZrZn 2 for po — 0.003 and initial values of 7 and 71 
which differ from 7 = 2"f± =0.2 only by numbers of order 
10~ 3 — 10~ 4 [26(. The obtained values for the coordinates 
of the maximum are: T m = 0.301 K at P m = 6.915 kbar. 
This result confirms the mentioned sensitivity of the lo- 
cation of the maximum T m towards slight variations of 
the material parameters. Experimental investigations of 
this low temperature-low pressure region with higher ac- 
curacy may help to locate this maximum with better pre- 
cision. 

We show by Fig. 3 the high- pressure part of same phase 
diagram . As far as the experimental phase diagram [l3j 



11 



has restricted accuracy in this range of temperatures and 
pressures, it is important to draw the attention to the 
first order phase transitions, shown by solid lines BC 
(FM-FS) and AC (N-FS) as this interesting topology of 
the phase diagram of ZrZn2 in the high-pressure domain 
[Pb < P < Pa) is not seen in the experiments. In fact 
the line AC is quite flat but not straight as the line BC. 
We hope that the detailed theoretical results for those 
parts of the phase diagram, where experimental data are 
lacking, may stimulate the further experiments in those 
ranges of pressure and temperature. The theoretical di- 
agram, shown in Fig. 1, looks almost the same as the 
experimental one as far as there the details are smeared. 
Note, that the theoretical diagram is a direct result of 
calculations with the free energy ([2|) and the proposed 
P-dependence of its parameters, without any approxima- 
tions and simplifying assumptions. 

The theoretical phase diagram reflects well the main 
features of the experimental behavior (TT{ , including the 
claimed change in the order of the FM-FS phase tran- 
sition at relatively high P. Within the present model 
the N-FM transition is of second order up to Pc ~ P c - 
Experiments give that the first order N-FM transition 
continues to much lower P values, and if this result is 
reliable, the theory can be easily modified to include this 
by a change of sign of bf. Then a new tricritical point 
will appear, located at some Pt r < Pc on the N-FM tran- 
sition line. Since Tc > 0, a direct N-FS phase transition 
of first order is predicted in accord with conclusions from 
de Haas-van Alphen experiments [32j an d some theoret- 
ical studies [HI . Such transition may not occur in other 
cases where Tc — 0. In spin-fluctuation theory(n = 2) 
the diagram topology remains the same but points B and 
C are slightly shifted to higher P, typically by about 
0.01 - 0.001 kbar. 



C. Diagram of UGe2 

The experimental data for UGe2 give Tfo = 52 K, 
P 9 = 1.6 GPa (= 16 kbarl, = 0.75 K, P m Pa 1.15 GPa, 

0, and P Qc pa 1.05 GPa HHEH- We US6 agahl 71 = 1 in 
the pressure dependence of Tf(P) From the above 
cited values for T m and P$ c we obtain 7 pa 0.0984 and 
71 ps 0.1678. The temperature Tc ~ 0.1 K corresponds 
to k — 4. 

In Fig. 4 we show the calculated with these initial pa- 
rameters T — P diagram of UGe2 . We have assumed that 
T s = in compliance with the arguments in Sec. III. A. 
The location of points A, B, C on the phase diagram is 
given by the following values: Ta = K, Pa = 1.723 
GPa, T B = 0.481 K, P B = 1.563 GPa, T c = 0.301 K, 
and Pc = 1.591 GPa. Figs. 5 and 6 show the details 
of low-temperature and the high-pressure parts of this 
phase diagram, respectively. 

Generally the main features of experimental curves 
are properly modelled by the theoretical T-P diagram 
of UGe2, although there is some discrepancy in the nu- 



merical values for P m corresponding to the maximum. 
Theoretically it is calculated as ~ 1.44 GPa in Fig. 4 
which is about 0.3 GPa higher than found experimen- 
tally 0, @. If the experimental measurement are pre- 
cise and give correct results for P m , this difference may 
be attributed to the so-called (T x ) metamagnetic phase 
transition in UGe2. It manifests in an abrupt change of 
the magnetization value in the vicinity of P m . The pos- 
sibility of metamagnetic transition is not considered in 
our model, but we may suppose that it will change the 
shape of Tps(P) making more symmetric the slope of 
curve around P m and lowering it together with Pb and 
Pc- We have tried to obtain a lower P m value, keeping 
T m unchanged within the model @, but this leads to a 
change of P c q to a value that disagrees with experiment. 
In spin-fluctuation approach (n = 2) the multi-critical 
points are located at slightly higher P (by about 0.01 
GPa), as for ZrZn2. Therefore, the results from spin- 
fluctuation theory are slightly worse than the results ob- 
tained from the usual linear approximation (n = 1) for 
the parameter t. 

In principle, the location of P m both for ZrZn2 and 
UGe2 may be tuned better within the present model |T]), 
if we take into account the anisotropy effects. Previous 
calculations [22j show that the inclusion, especially of 
spin-triplet Cooper pair anisotropy (the u s -term in the 
Eq. |1|), changes the values of coordinates of points A, 
B, C and max. The problem is that the anisotropy in the 
present status of knowledge may serve only as a fitting 
parameter, because no proper information from micro- 
scopic theories is available both for its sign and magni- 
tude. But such fitting from general considerations can 
obscure other effects, as the mentioned above metamag- 
netic transition in UGe2 , and even lead to some wrong in- 
terpretations when comparing the experimental and the- 
oretical curves. 



D. Two types of ferromagnetic superconductors 
with spin-triplet electron pairing 

The estimates for UGe2 give that 71 k pa 1.9, so the 
condition for Tps(P) to have a maximum found from 
Eq. ([12")) is satisfied. As we have discussed for ZrZii2, 
the location of this maximum can be hard to fix accu- 
rately in experiments. However, P c q can be more easily 
distinguished, as in the UGe 2 case. Then we have a well- 
established quantum (zero-temperature) phase transition 
of second order, i.e., a quantum critical point at some 
critical pressure Pq c > 0. As shown in Sec. III.C, un- 
der special conditions the quantum critical points could 
be two: at the lower critical pressure Pq c < P m and the 
upper critical pressure Pq c < P m . This type of behav- 
ior in systems with T s = (as UGe2) occurs when the 
criterion (fl"8|) is satisfied. Such systems (which we label 
as U-type) are essentially different from those such as 
ZrZn2 where 71 < 7 and hence Tps(0) > 0. In this lat- 
ter case (Zr-type compounds), a maximum T m > may 
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sometimes occur, as discussed earlier. We note that the 
ratio 7/71 reflects a balance effect of two interactions be- 
tween xj}-M. When the trigger interaction, represented 
by 7, prevails, the Zr-type behavior is found where su- 
perconductivity exists at P = 0. The same ratio can be 
expressed as jo/dMo, which emphasizes that the ground 
state value of the magnetization at P = is also rel- 
evant. Alternatively, one may refer to these two basic 
types of spin-triplet FSs as "type I" (for example, for the 
"Zr-type compounds"), and "type II" - for the U-type 
FS compounds. 

As we see from this classification, the two types of spin- 
triplet ferromagnetic superconductors have quite differ- 
ent phase diagram topologies although some fragments 
have common features. The same classification can in- 
clude systems with T s ^ but in this case one should 
use the more general criterion (|24p. 



E. Other compounds 

In URhGe, 7/(0) ~ 9.5 K and T FS (0) = 0.25 K and, 
therefore, as in ZrZn2,here the spin-triplet superconduc- 
tivity appears at ambient pressure deeply in the ferro- 
magnetic phase domain [f| @, H| . Although some similar 
structural and magnetic features with UGe2 the results 
111 Ref. n of measurements under high pressure show 
that, unlike the behavior of ZrZn2 and UGe2, the ferro- 
magnetic phase transition temperature Tp(P) ~ Tf(P) 
has a slow linear increase up to 140 kbar without any 
experimental indications that the N-FM transition line 
may change its behavior at higher pressures and show a 
negative slope in direction of low temperature up to a 
quantum critical point Tp — at some critical pressure 
P c . Such behavior of the generic ferromagnetic phase 
transition temperature cannot be explained by our ini- 
tial assumption for the function Tf(P) which was in- 
tended to explain phase diagrams where the ferromag- 
netic order is depressed by the pressure and vanishes at 
T = for some critical pressure P c . The Tps(P) line of 
URhGe shows a clear monotonic negative slope to T = 
at pressures above 15 kbar and the extrapolation [8| of 
the experimental curve Tps(P) tends a quantum critical 
point T F s(Po C ) = at P 0c - 25 - 30 kbar. Within the 
framework of the phenomenological theory ©, this T—P 
phase diagram can be explained after a modification of 
X/(P)-dependence is made, and by introducing a conve- 
nient nontrivial pressure dependence of the interaction 
parameter 7. Such modifications of the present theory 
are possible and follow from important physical require- 
ments related with the behavior of the /-band electrons 
in URhGe. Unlike UGe2, where the pressure increases 
the hybridization of the 5/ electrons with band states 
leading to a suppression of the spontaneous magnetic mo- 
ment M, in URhGe this effects is followed by a stronger 
effect of enhancement of the exchange coupling due to the 
same hybridization, and this effect leads to the slow but 
stable linear increase of the function T F (P) @|. These 



effects should be taken into account in the modelling of 
the pressure dependence of the parameters of the theory 
©when applied to URhGe. 

Another ambient pressure FS phase has been observed 
in experiments on UCoGe [9j. Here the experimentally 
derived slopes of the functions Tp(P) and Tps(P) at rel- 
atively small pressures are opposite compared to those 
for URhGe and, hence, the T — P phase diagram of this 
compound can be treated within the present theoretical 
scheme without substantial modifications. 

Like in UGe2, the FS phase in Ulr [12j is embedded in 
the high-pressure/low-temperature part of the ferromag- 
netic phase domain near the critical pressure P c which 
means that Ulr is certainly a U-type compound (see Sec. 
IV. D). In UGe2 there is one meta- magnetic phase transi- 
tion between two ferromagnetic phases (FM1 and FM2), 
in Ulr there are three ferromagnetic phases and the FS 
phase is located in the low-T/high-P domain of the third 
of them - the phase FM3. There are two meta-magnetic- 
like phase transitions: FM1-FM2 transition which is fol- 
lowed by a drastic decrease of the spontaneous magneti- 
zation when the lower-pressure phase FM1 transforms to 
FM2, and a peak of the ac susceptibility but lack of ob- 
servable jump of the magnetization at the second (higher 
pressure) "meta-magnetic" phase transition from FM2 to 
FM3. Unlike the picture for UGe2, in Ulr both transi- 
tions, FM1-FM2 and FM2-FM3, are far from the max- 
imum T m (P TO ), so in this case one can hardly speculate 
that the max is produced by the nearby jump of magne- 
tization. Ulr seems to be a U-type spin-triplet ferromag- 
netic superconductor. 



V. FINAL REMARKS 

Finally, even in its simplified form, this theory has been 
shown to be capable of accounting for a wide variety of 
experimental behavior. A natural extension to the theory 
is to add a M 6 term which provides a formalism to inves- 
tigate possible meta-magnetic phase transitions [37j and 
extend some first order phase transition lines. Another 
modification of this theory, with regard to applications to 
other compounds, is to include a P dependence for some 
of the other GL parameters. The fluctuation and quan- 
tum correlation effects can be considered by the respec- 
tive field-theoretical action of the system, where the order 
parameters ij) and M are not uniform but rather space 
and time dependent. The vortex (spatially non- uniform) 
phase due to the spontaneous magnetization M is an- 
other phenomenon which can be investigated by a gener- 
alization of the theory by considering non-uniform order 
parameters fields ip and M (see, e.g., Ref. [38, 39]). Note, 
that such theoretical treatments are quite complex and 
require a number of approximations. As already noted in 
this paper the magnetic fluctuations stimulate first order 
phase transitions for both finite and zero phase transition 
temperatures. 
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